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ABSTRACT 

I explore the hypothesis that many of the unusual aspects of the apparent stellar bar of the Large 
Magellanic Cloud are the result of viewing a triaxial stellar bulge that is embedded in a highly 
obscuring thin disk. Specifically, this hypothesis explains the observed off-center position of the bar 
within the disk, the differing apparent distances of the bar and disk, the near alignment of the bar's 
major axis position angle and the disk's line-of-nodes, and the asymmetric appearance of the bar 
itself. Indirectly, it may also play a role in explaining the microlensing rate toward the LMC and the 
recently observed large velocity dispersion of RR Lyrae stars. 

Subject headings: Magellanic Clouds — galaxies: photometry — galaxies: stellar content 



1. INTRODUCTION 

The off-center bar of the Large Magellanic Cloud 
LMC) has helped define a class of irregular galaxies 
de Vaucouleurs Freemanlll972D and become a prime 
motivation for the study of off-center, spiral wa ve driv- 
ing m odels ijColin fc Athanassoula 1 989; Gardiner et al.l 
[ilia)- Despite the obvious "barred" appearance of the 
LMC (Figure^), there are surprisingly few direct con- 
straints on the t hree dimensional stru cture of the appar- 
ent feature fsee lZhao fc Evan3l2000|) . Is the bar in the 
LMC real or an illusion? 

Several features of the LMC bar are difficult to recon- 
cile with the simplest version of barred disk galaxies: 

1) The bar is off-center relative to the underlying disk. At 
large radius (i? > 5 kpc), where there is no contamina- 
tion of the disk isophotes by other components and where 
the deprojected position angle, ellipticity, and centroid 
remain approximately constant, the disk center is mea- 
sured to lie a pproximately at a = 5''29'"and — 69°30' 
(see Fig 4 of Ivan der Marell 1200]]) . At small to inter- 
mediate radii (0.5 kpc < R < 2.5 kpc), where the bar 
dominates the isophotes and the deprojected position 
angle and centroid remain approximately constant, the 
centroid is shifted from that of the ou ter disk by ~ 
0.5 a nd 0.3 kpc in x and y, respectively ijvan der Marell 
1200 This type of misalignment has been defined to 
be a common characteristic of barred Magellanic irregu- 
lars (|de Vaucouleurs fc Freemanlll972>) and dynamically 
important in defining the weak spiral structure of these 
systems, but its cause has not yet been conclusively iden- 
tified (Gardiner et al. 1998). 

2) The meas ured distance to the b ar is different than that 
to the disk. iNikolaev et alJ l)2004|) conclude that the bar 
is "elevated" above the disk plane, located closer to us 
by ~ 0.5 kpc. 

3) The bar's major axis is nearly parallel to the disk's 
line-of-nodes ( 1 25° ± 5° and 122.5° ± 8.3°, respectively; 
Ivan der Marell 1)2001(1 ). The foreshortening of any bar 
whose major axis position angle is not oriented along the 
line-of-nodes will make the alignment appear closer than 
it truly is, but most bars in disk galaxies are sufficiently 
elongated that they are typically not aligned in projec- 
tion. 

4) The bar is asymmetric. The stellar density is ellip- 



soidal on the southwestern side and flat along the north- 
eastern edge. The flat edge lies nearly parallel to the 
disk's line-of-nodes (Figure QJ)). 

The hypothesis I explore is that the LMC's barred ap- 
pearance and the features listed above are the result of 
viewing a bulge population that is obscured by an opti- 
cally thick inner disk. 

2. A SIMPLE MODEL 
2.1. The Observational Situation 

The data come from the Magellanic Clouds Photomet- 
ric Survey, which is a J7, -B, V, and / survey of the central 
8.5° X 7.5° of the Large Magellanic Cloud and 4.5° x 4° of 
the Small Magellanic Cloud . Both catalogs are published 
(|Zaritskv et alj[200ll200l . 

I plot in Figure^the stellar density on the sky of stars 
with 16 < V < 20. The high-contrast panels (Figures 
nj3,d), show how the distribution of the central popula- 
tion appears to be sharply truncated toward the north- 
east and rounded toward the southwest. The line-of- 
nodes and center derived by van der Marel (2001) are 
plotted for reference in Figure Some rectilinear ir- 
regularities in the density distribution that are present in 
Figure m arise from different completeness limits among 
drift scan images in this region of high stellar density. 
Although these irregularities do not affect the qualita- 
tive conclusions presented here, they do preclude detailed 
model fitting to these data prior to obtaining the results 
from extensive artificial star tests. I do, however, provide 
xi values for the various models presented for intercom- 
parison, assuming only Poisson noise in the star counts. 
These values are calculated within a rectangular region 
than incorporates the visible bar in Figure and the 
model is normalized to have the same total number of 
stars within this region. Because of the systematics de- 
scribed above, and the lack of a modeled disk component, 
xt is expected to be 1 even for the correct model. 

2.2. Exploring the Hypothesis 

I begin wi th a m odel defined by the disk geometry 
measured bv lvan de r Marel & Cioni (2001). The choice 
of geometry, in particular the position angle of the line- 
of-nodes (122° in that study), will be explored below. 
Within the adopted geometry, the bulge has two of its 



2 



N 



E 




Fig. 1. — The observed and modeled stellar density distribution in the LMC. Panel (a) is a low contrast map of the stellar density for 
stars with 16 < V < 20. Panel (b) is a high contrast map of the same stellar density distribution. There is a central population that 
is sharply truncated toward the northeast, rounded toward the southwest, and roughly aligned with the disk's line-of-nodes. Panel (c) 
shows the distribution of bulge particles in the model described in the text. Panel (d ) is the same high contrast map shown in (b) with an 
overlay indicating the geometry of the disk and bulge system. The line-of-nodes ffrom iran der Marel Hi Cioni 1 1200 ID ') is represented by the 
straight line that bisects the rectangle, and the position of LMC center as measured from the outer isophotes is marked with a tick mark 
perpendicular to the line-of-nodes that is slightly off-center relative to the ellipsoid that delineates the bulge appearance in the absence of 
the optically thick disk. 



principal axes in the disk plane (and hence the third 
perpendicular to the disk plane), but has an otherwise 
unconstrained position angle and axial ratio. Its cen- 
ter is assumed to be coincident with that of the disk, 
because the focus here is to determine whether the ob- 
servations can be explained without invoking an offset 
component. I assume that all of the obscuring dust 
is in the disk plane and has a much smaller vertical 
scaleheight than any dynamically hot population (see 
[Zaritsky et al. 2004, for a confirmation of this assump- 
tion for the LMC as a whole). The volume density of 
bulge stars is assumed to be constant for 0.001 < r < 0.3, 
where r = y/ (x/a)^ + [y/by + {z/cY, to follow an r^^ 
density profile for 0.3 < r < 4, and to otherwise be zero. 
The scales, a, h, and c, are free parameters. The constant 
density core is set to avoid a divergence at small radii and 
to mimic the observed shallow density gradient in the in- 
ner kpc. Because of crowding in the MCPS images, the 
true profile may be significantly steeper than that sug- 
gested by Figure n and hence, the flat central density of 
the models may not reflect reality even though it is ulti- 
mately successful in reproducing the observations. I draw 
stars from a power law luminosity function, N{ra) oc m? , 
place them in projection according to the viewing geom- 



etry of the LMC, apply extinction (extinction prescrip- 
tion is discussed below) , and accept only those stars that 
satisfy the observational criteria. Distance differences 
are accounted for in determining the apparent magni- 
tude and no disk component is included in the generated 
stellar density maps. In summary, the parameters of the 
model that vary are the three scaling parameters, a, 6, 
and c, the core radius in scale-free units (although it was 
fixed for all models once a satisfactory value was found), 
and the bulge's position angle. 

After exploring various parameter combinations, I 
reached a few principal conclusions. First, the extinc- 
tion in the midplane, at least in this central portion 
of the LMC disk, must be significantly larger than es- 
timated by Zaritskv et al. ( 2004) to sufficiently obscure 
the rear portion of the bulge. To increase the extinction 
in the modeling, I define a minimum visual extinction 
value, Amin- At any location for w hich the extinction 
value from the iZaritskv et alJ l)2004|) extinction map is 
less than this minimum value, the extinction is set to 
Amin- An extreme value of Amin = 10 was used in Fig- 
ure ^ to mimic complete extinction, but a strong ex- 
tinction edge is seen toward the northeast side of the 
bulge for Amin > 2. These large extinction values do not 
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Fig. 2. — The effect of changing the y and z bulge scales. Panel 
(a) shows a model with a = 2.5, b = 1.0, and c = 2.5. Panel (b) 
shows a model with a = 2.5, b = 0.5, and z = 1.5. For comparison, 
the model in Figure [Tl has a = 2.5, fe = 0. 5 and c = 2.5. All of 
these models adopt the Ivan der Marel Ik C ioni 12001) geometry. 
The image is cropped to only show the inner region of the LMC. 



contradict the lZaritskv et all l|2004j) results because such 
highly obscured stars would not be in their sample. Sec- 
ond, the extension of the bulge must be relatively short 
along the y-axis (the disk plane, short axis of the bulge) 
because otherwise the sharp northeastern edge created 
by the disk bisection of the bulge becomes increasingly 
blurred (see Figure for a model with double the y- 
scale of that shown in Figure H^). Third, the extension 
of the bulge along the z-axis (out of the disk plane) must 
be relatively large to produce the noticeable rounding of 
the southwestern edge (Figure EJd). 

The parameters for the model shown in Figure^ are 
a = 2.5 kpc, b — 0.5 kpc, and c = 2.5 kpc, which makes 
this bulge an oblate spheroid that protrudes from the 
disk plane (this model has xl = 9-9j while the same 
geometric model without any midplane extinction has 
xl = 46.9). The physical scale of the constant density 
central region along each axis is given by the core bound- 
ary, defined previously to be r = 0.3, times either a, b, 
or c. The extreme nature of this bulge can be mitigated 
somewhat, at the expense of a moderate worsening of 
the fit. For example, a bulge with a — 2.5, b — 1.0 and 
c = 2.5 still reproduces the principal qualitative struc- 
ture but does not produce as sharp an edge toward the 
northeast (Figure [5^). This model has a slightly lower 
xi, 9.2, than the previous model presumably because it is 
more extended and so partially accounts for the missing 
disk stars in the model. The extent along the z axis can- 
not be decreased significantly without losing the rounded 
appearance of the bar (see Figure [2b for a model with 
a = 2.5, b = 0.5, and c = 1.5; xl = 76.8). Any detailed 
exploration of the parameter space would require a disk 
model in which the bulge is embedded, and higher reso- 
lution imaging to produced stellar density distributions 
that are less severely affected by crowding. 

Finally, I explore the effect of changing the adopted 
position angle of the disk's line-of-nodes. The ini- 
tial angl e in the models was that measured by 
Ivan der Marel fc Cioni I lj200lD , but another recent mea- 
surement produced an angle that is ^ 30° larger 
ijNikolaev et al.ll200^ . Surpr isingly, this difference does 
not affect the model's ability to recreate the bar appear- 



ance. By changing the bulge position angle from being 
10° off the line-of-nodes to 45°, I reproduce the observed 
appearance (Figure |31 xl = 10.1, a = 3.0, & = 0.5, and 
c = 3.0). This new geometry helps explain one feature 
that was not well accounted for in the original models. 
The circles in Figure O represent the center of the LMC 
as determined from the outer disk isophotes. The ob- 
served midpoint of the "bar" along the line-of-nodes is 
displaced by 15 arcmin toward the northwest, while the 
stellar density weighted centroid of the model presented 
here is displaced by 23 arcmin in the same northwesterly 
direction. In the simulations, at least, this offset from 
the center is due to the longer line of sight through the 
bulge to the disk on one side of the LMC center than 
on the other. The uncertainties in determining a pre- 
cise LMC center and in the stellar density distribution 
of the MCPS in this highly crowded region preclude one 
from inverting the model fits to reach a conclusion on the 
correct position angle of the line-of-nodes. Nevertheless, 
is is possible to reproduce the qualitative nature of the 
observations over the range of observed values for the po- 
sition angle of the disk's line-of-nodes. The relative ease 
with which the general model was able to fit the data 
despite a large change in the line-of-nodes demonstrates 
that its ability to produce populations that mimic the 
observations is not predicated on a special coincidence of 
paramaters. 

3. DISCUSSION AND SUMMARY 

Although this modeling exercise demonstrates that the 
LMC might have a stellar bulge that when viewed in pro- 
jection across an optically thick disk produces a feature 
that appears to be an off-center, levitating bar, it does 
not prove that this is indeed the case. The principal 
differences between a bar and a bulge in the LMC lie 
in the line-of-sight extent and the vertical (out of the 
disk plane) velocity dispersion. There are tantalizing ob- 
servations that suggest that these might both be large, 
and hence support the bulge hypothesis. First, aside 
from the different distance measurement to the disk and 
bar (Nikolacv et al. 2004), a protruding stellar popula- 
tion would help reproduce the lar ger-than-expected mi - 
crolensin g rates toward the LM C ijZhao Sz Evanjl2000() . 
Although lZhao fc Evanal 1)20001) speculated that the lens- 
ing arose from an unrelaxed stellar population that 
in projection appeared to be a bar, the hypothesized 
bulge will produce an analogous lensing enhancement 
and the strong disk-plane extinction would address the 
spatial asymmetry of the discovered microlensing events 
ijMancini et alJl2'004tl . Whether this hypothesis survives 
further, more quantitative, scrutiny depends on the mass 
and radial profile of the bulge, both of which are un- 
known j_bu^_evaJimtmg^^ similar to those carried 
out bv lMancini et alJ 120041) for this newly proposed ge- 
ometry is the next step. 'Alv§^ ()2004|) has fit King models 
to the distribution of RR Lyrae stars, which presum- 
ably trace an old population, and finds a core radius of 
1.42 ± 0.12 kpc (a rough match to the core radii of 0.7 
to 1.0 kpc imposed in these models), and a steep decline 
(exponential scale length of 1.47 ± 0.08 kpc). This mea- 
surement might be providing a description of the radial 
distribution of the bulge population. Second, there is ev- 
idence for^^^inBnmtkallyho^ in the 
LMC ijMinniti et al.ll2003l: iBorrisova et al.ll2004() . Those 
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Fig. 3. — Comparison of model (a) an d data (b) if the dis k line-of-nodes is at 150°. Filled circles denote the center of the LMC as 
determined from the outer disk isophotes Ivan der "W^eH^^. The image is cropped to only show the inner region of the LMC. 



authors suggest that this component is associated with 
the stellar halo, but because they observe the high veloc- 
ity sta rs at small proje cted radii (for example, R < 1.5° 
for theH orrisova et al.l l)2004|) study) , rather than at radii 
where only LMC "halo" would exist, these stars may 
instead be members of a dynamically hot bulge com- 
ponent. The large velocity dispersion ('^ 53 km/sec; 
iBorrisova et al. 2004) implies a significantly larger ver- 
tical scaleheight than that of the disk for this population 
(the kinematically hottest disk component has a velocity 
dispersion of ^ 20 km/s ec; Meathcringham et al. 1983 
Ivan der Marel et al.ll2002i) . 

Now I revisit each of the "anomalous" observations dis- 
cussed in §1 in the context of the bulge plus obscuring 
disk hypothesis: 

1) The apparent bar is off-center because we only see 
half the bulge, and that half at an odd angle. The gen- 
eral geometrical features of the appar ent bar are repro- 
d uced well with either th e fvan der Mar el & Cioni (200lJ) 
or lNikolaev et alJ l|2004j) disk geometry (Figuresnand|31) • 
The latter does slightly better at fully reproducing the 
apparent misalignment of disk and bar centroids. 

2) The bar and disk appear to be at different distances 
because the bar stellar sample is biased in favor of stars in 
front of the disk. The apparent average displacements in 
distance of the entire bar population relative to the disk 
center for the models presented in Figures^ and|3ti' are 
1 and 0.6 kpc, respectively. T hese are both larger th an 
the measurement presented bv iNikolaev et al"l l)2004f) of 



^ 0.5 kpc, but there are various simple ways in which 
the models could reproduce the published value. For 
example, if the mid-plane extinction is reduced slightly 
then some stars behind the disk will enter the sample. 

3) The bar major axis is nearly aligned with the disk 
line-of-nodes because the obscuring sheet cuts through 
the bulge along the line-of-nodes. As such, the align- 
ment of the projected major axis with the line-of-nodes is 
less problematic in this obscure d bulge model than in an 
unobscured bar model. For the Ivan der Marel fc Cionil 
l)2001|) geometry, the bulge major axis in our model 
is within ^ 1 0° o f the line-of-nodes, but for the 
INikolaev et al.l lj2004jl geometry the difference in posi- 
tion angles is ~ 45°. Within the range of observationally 
allowed line-of-node angles, the model does not require 
precise alignment between the line-of-nodes and bulge 
position angle to reproduce the observations. 

4) The apparent bar asymmetry arises because the north- 
eastern side is artificially truncated by the obscuring disk 
(Figure HI. 

Despite these successes, the bulge hypothesis is not 
without challenges. First, the shape of the derived 
bulge, with a long or intermediate axis out of the plane, 
is unusual. This geometry is driven by having to re- 
produce the extent of the bar in the southwestern di- 
rection. If the disk plane is actually more highly in- 
clined to the line of sight than in the models, then 
the bulge would not need to extend as far out of the 
disk to project as far on the sky. Because the inch- 
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nation is well determined for the outer disk, such an 
adjustment in inclination would require the LMC disk 
to be warped. Second, the bar feature does not ap- 
pear grossly differen t in infrared stellar density images 
ijvan der M arcl 2001"). For extinction to strongly affect 
the 2/xm isopohotcs would require an optical extinctio n 
> 10 {AkIAv ~ O.lr lCardelh. Clayton, fc Mathi^lTosl . 
Such a large, rather uniform extinction appears to be un- 
likely. Again, a more highly inclined disk would help the 
situation because the path length through the disk would 
be longer and the extinction wou ld be larger. Lastly , 
existing reddening measurements l|Nikolaev et al. 1 120041: 
iZaritskv et alJ2004D do not find a significant tail of mod- 
erately extincted (1 < Ay < 2) stars. Unfortunately, the 
iNikolaev et al.l l|2004|) Cepheid sample is sparse in the 
region of interest immediately northeast of the bar (con- 
sisting of only a few tens of stars), the IZaritskv et aTl 
l)2004|) photometry is strongly affected by crowding in 
the inner LMC, and the extinction may be sufficient to 
eliminate the highly extincted stars from these samples 
in the first place. 

How can we make further progress in determining 
whether the obscured bulge hypothesis is correct? A 
complete kinematic survey of the LMC from which one 



could determine whether the dynamically hot compo- 
nent has a similar spatial distribution as the projected 
bulge and whether, or not, there are kinematic signa- 
tures within the disk of a bar pote ntial could provide 
a definitive test. iZhao et'aP l|2003^ present an exten- 
sive kinematic survey, but for this purpose their velocity 
resolution is too coarse (15 km sec~^) and the bar re- 
gion is undersampled. Infrared surveys should eventually 
penetrate through the dust and provide evidence for or 
against a hidden bulge population. And so, on the basis 
of the available data, I can only suggest that the conun- 
drums posed by the LMC bar may be resolved by the 
presence of a triaxial bulge and an optically thick inner 
disk, and look forward to learning whether this is indeed 
the case. 
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